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Abstract. It is shown that amplification of the magnetic field in supernova remnants (SNRs) occurs in all six objects 
where morphological measurements are presently available in the hard X-ray continuum at several keV. For the 
three archetypical objects (SN 1006, Gas A and Tycho's SNR) to which nonlinear time-dependent acceleration 
theory has been successfully applied up to now, the global theoretical and the local observational field strengths 
agree very well, suggesting in addition that all young SNRs exhibit the amplification effect as a result of very 
efficient acceleration of nuclear cosmic rays (CRs) at the outer shock. Since this appears to be empirically the 
case, we may reverse the argument and consider field amplification as a measure of nuclear CR acceleration and it 
has indeed been argued that acceleration in the amplified fields allows the GR spectrum from SNRs to reach the 
knee in the spectrum or, in special objects, even beyond. The above results are furthermore used to investigate the 
time evolution of field amplification in young SNRs. Although the uncertainties in the data do not allow precise 
conclusions regarding this point, they rather clearly show that the ratio of the magnetic field energy density and 
the kinetic energy density of gas fiow into the shock is of the order of a few percent if the shock speed is high 
enough Vb > 10^ km/s, and this ratio remains nearly constant during the SNR evolution. The escape of the highest 
energy nuclear particles from their sources becomes progressively important with age, reducing also the cutoff in 
the TT^-decay gamma-ray emission spectrum with time after the end of the sweep-up phase. Simultaneously the 
leptonic gamma-ray channels will gain in relative importance with increasing age of the sources. 

Key words, cosmic rays, shock acceleration, nonthermal emission, supernova remnants; individual: Tycho's super- 
nova remnant 



1. Introduction 

Recent observations with the Chandra and XMM-Newton 
X-ray space telescopes have confirmed earlier detec- 
tions of nonthermal continuum emission in hard X-rays 
from young shell- type supernova remnants (SNRs). With 
Chandra it became even possible to resolve spatial scales 
down to the arcsec extension of individual dynamical 
structures like shocks. Results of this type have been 
pubhshed for Tycho's SNR (G120.1-bl.4) (Hwang et al. 
2002; Bamba et al. 2003b), RCW 86 (G315.4-2.3) (Rho 
et al. 2002; Bamba et al. 2003b), Cas A (Gill. 7-2.1) 
(Vink & Laming 2003; Bamba et al. 2003b), SN 1006 
(G327.6-fl4.6) (Long et al. 2003; Bamba et al. 2003a; 
2003b; 2004), and Kepler's SNR (G4.5-h6.8) (Bamba et 
al. 2003b); (see also Bamba 2004). Small-scale filamentary 
structures in SNR RX J1713. 7-3946 (G347.3-0.5) have 
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been reported by Uchiyama et al. (2003) and Lazendic et 
al. (2004). 

These results can be seen from three different view- 
points: First of all, assuming the filamentary structures to 
represent synchrotron emission near shocks, their width 
can be considered as the synchrotron coBoling length to 
derive the effective downstream magnetic field strength 
-Bd- Secondly, we can compare this with the value for 
the interior magnetic field strength that makes a com- 
plete dynamical model of the SNR consistent with the 
spatially integrated synchrotron spectrum, from the ra- 
dio range to hard X-rays. The third viewpoint is from the 
plasma physics of field amplification in an ionized gas with 
a very strong nuclear cosmic ray (CR) gradient. 

Let us first turn to the observed X-ray morphology. 
In their work on Cas A, Vink & Laming (2003) have in- 
terpreted some of the filamentary hard X-ray structures 
as the result of strong synchrotron losses of the emitting 
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multi-TeV electrons in amplified magnetic fields down- 
stream of the outer accelerating SNR shock. Such an in- 
terpretation has been independently given by Berezhko et 
al. (2002) from a consistent description of the observed 
SN 1006 dynamics and of the properties of its overall non- 
thermal emission, and by Berezhko et al. (2003a) for the 
Chandra observations of SN 1006, leading to an effective 
downstream field strength of 120 /iG for the latter ob- 
ject. Subsequently Berezhko & Volk (2004a) revised the 
Vink & Laming result for Cas A by consistently taking 
account of all relevant dynamical and geometrical (pro- 
jection) factors. It was shown that for the highest energy 
electrons, usually responsible for the X-ray emission (like 
in the case of SN 1006 and Cas A), synchrotron losses be- 
come important already during the acceleration process. 
Together with the projection effect they imply a signifi- 
cantly shorter electron cooling time behind the shock, and 
thus an effective downstream magnetic field strength in 
Cas A that is stronger by a factor of five than estimated 
by Vink & Laming (2003). It reaches about 0.5 mC in the 
immediate downstream region of the outer shock. 

Such observational results gain their qualitative signif- 
icance through the fact that for SN 1006 (Berezhko et al. 
2002) and Cas A (Berezhko et al. 2003b) these effective 
magnetic fields and morphologies turned out to be exactly 
the same as predicted theoretically from nonlinear acceler- 
ation theory. The theory uses gas dynamics and nonlinear 
acceleration theory to calculate a kinetic model of young, 
evolving SNRs which includes the energetic electron com- 
ponent together with the nuclear CR component. Since 
this theory does not explicitly contain equations for the 
evolution of the electromagnetic fields or for the injection 
of low-energy particles into the acceleration process, these 
processes are parametrized by an effective magnetic field 
strength and an injection strength for the nuclear par- 
ticles. The theory then works as follows: the comparison of 
the model with the observed spectral features of the spa- 
tially integrated synchrotron emission from the remnant 
yields a selfconsistent value for Bd (see also Reynolds & 
Ellison 1992) as well as for the total CR pressure Pc- The 
amplitude of the observed synchrotron emission limits in 
addition the pressure of the energetic electrons to 1 per- 
cent or less of the total CR pressure Pc in the remnant. 

It is important to note that the magnetic field B'^ de- 
rived from the fit of the calculated overall synchrotron fiux 
to the observed flux and the fleld B'J^ that results from 
the observed X-ray brightness profile correspond to basi- 
cally different parts of the SNR. The field belongs to 
the thin region just behind the shock whereas P^, which 
mainly comes from the radio data, characterizes the av- 
erage effective field within a much wider region, occu- 
pied by the GeV-electrons. Generally speaking, B[^ and 
Pj can therefore also be basically different. The equality 
P^ = Pj is expected if during the SNR evolution the ran- 
dom amplified field is generated near the shock front at 
the level B'^^ oc V^^ (Berezhko et al. 2003a; Berezhko & 
Volk 2004b). Therefore the equality of the experimental 



values P^ and P^' confirms that the field is indeed ampli- 
fied (generated) to the above level. 

It should perhaps be noted here that there is a good 
semi-quantitative understanding of the injection rate of 
nuclear particles into the diffusive shock acceleration pro- 
cess in terms of the escape of suprathermal downstream 
particles. However, due to the incompletely known relax- 
ation state of the downstream ions, the escape rate is the- 
oretically known only up to factors of order unity (Volk et 
al. 2003), and therefore it is necessary to fix its value em- 
pirically in a quantitative model for the SNR as a whole. 
Overall, with the three parameters Pj, Pc, and electron- 
to-proton ratio determined empirically, it is possible to 
describe the overall observed electron synchrotron spec- 
trum from first principles. 

The large total CR pressure and the negligible contri- 
bution to it from electrons is a compelling argument for 
the existence and dominance of the nuclear CR compo- 
nent in SNRs which makes up almost the entire pressure 
Pc- The implied magnetic field strength must be the same 
as that determined independently from the observed X-ray 
filaments. This is a direct empirical test for acceleration 
theory. Even though such a detailed agreement has been 
demonstrated only in two cases up to now, it is hard to 
believe that this is a mere coincidence. 

Besides morphology and acceleration theory the third 
aspect of field amplification concerns plasma instability 
theory. In shocks with strong nuclear CR production the 
downstream CR pressure Pc is of the order of the ram 
pressure pV^ of the upstream gas flow in the shock frame 
(Drury & Volk 1981; Axford et al. 1982). Then the unsta- 
ble streaming of the shock accelerated nuclear particles 
(BeU 1978; Blandford & Ostriker 1978) leads to an ex- 
cessive resonant excitation of Alfvcn waves whose energy 
density exceeds that of the mean magnetic field by far 
(McKenzie & Volk 1982), making an increased and max- 
imally disordered magnetic field plausible (Volk 1984). 
Much more recently, analytical studies and numerical sim- 
ulations of this instability were performed by Lucek & Bell 
(2000) and Bell & Lucek (2001) which indicated strong 
field amplification for Pc ~ pV^, while the spatial diffusion 
for energetic particles in this disordered field reached the 
limiting level of Bohm diffusion. This implies that the scat- 
tering mean free path approaches the gyro radius in the 
amplified field. Finally, Bell (2004) extended an early in- 
vestigation by Achterberg (1983) of the nonresonant right- 
hand polarized low frequency MHD mode, propagating 
parallel to the average magnetic field in the presence of a 
nuclear CR current. In the case Pc ^ pV^ he found strong 
linear growth with a growth rate in excess of the oscilla- 
tion frequency, i.e. a purely growing mode. This growth is 
even more rapid than that of the resonant Alfven waves, 
usually considered until now. Nonlinear quasi-MHD simu- 
lations on the basis of this instability showed strong field 
amplification for the energetic particle streaming that one 
has to expect in the strong accelerating shocks of young 
shell-type SNRs. 
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We assume the observed filamentary X-ray structures 
to be manifestations of sucli tlicoretical plasma processes 
and shall therefore consider them as an integral aspect 
of diffusive shock acceleration in violent events such as 
supernova explosions. 

The aim of our paper is to discuss the objects where 
field amplification possibly occurs. Therefore we shall con- 
sider all those shell-type SNRs where filamentary shock 
structures in hard X-rays have been published up to now. 
These SNRs are typically quite young, with ages not much 
in excess of a sweep-up time for a uniform circumstellar 
medium. The explanation for this selection effect is as fol- 
lows: the cutoff frequency of synchrotron spectrum pro- 
duced by the shock accelerated electrons is proportional 
to V^^ and diminishes after sweep-up as a result of the 
slowing down of the shock, reducing the X-ray emission 
at a fixed observation energy energy, like 4 keV, sharply 
with time (Berezhko & Volk 2004b). This is independent 
of the fact that the total CR energy, the normalized CR 
pressure Pc/ {pV^), and the hadronic gamma-ray luminos- 
ity remain quite high for much longer times in the Sedov 
phase (e.g. Berezhko & Volk 1997). Older SNRs are there- 
fore not expected to exhibit significant nonthermal hard 
X-ray emission. 

To search for field amplification in young SNRs only 
does therefore not imply a methodological restriction. On 
the contrary, if wc knew that all young SNRs showed field 
amplification we would know that all SNRs have a strong 
nuclear energetic particle population, as a result of the 
strong theoretical connection between field amplification 
and the dominant contribution of nuclear CRs to the over- 
all nonthermal pressure Pc- 

Besides Gas A and SN 1006 the available Chandra ob- 
servations concern Tycho's and Kepler's SNRs, RCW 86, 
and RX J1713.7-3946, as mentioned above. Tycho's SNR 
plays a special role, because there exists a detailed theo- 
retical model for the particle spectra and the synchrotron 
morphology (Volk et al. 2002). We shall use the projec- 
tion model for the brightness morphology developed by 
Berezhko & Volk (2004a) wherever possible to fit the data 
by adjusting the strength of the effective magnetic field 
i?d. Wc find that all the investigated sources have ampli- 
fied fields to a varying degree, and therefore we propose 
that all these sources are efficient nuclear CR accelerators. 
This is a necessary condition for the Galactic SNRs to con- 
stitute the source population of the Galactic CRs, at least 
up to the so-called knee of the energy spectrum. Regarding 
Tycho's SNR we shall find that it constitutes the third case 
where spectral and morphological synchrotron character- 
istics independently give the same effective field strength. 



2. Effective magnetic fields from X-ray data 

As shown by Berezhko & Volk (2004a), the width L « 
7^2 of the observable brightness profile J,^{e^,p) is always 
appreciably wider than the downstream width I2 in radius 



r of the three-dimensional emissivity q2{€^,r), simply for 
geometric reasons: 
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In this approximate relation, valid for I2 ^ O.lPg, Ra and 
p are the shock radius and the distance between the cen- 
ter of the remnant and the line of sight, respectively. At 
the same time I2 is the distance through which the very 
energetic electrons stream away from the shock during a 
synchrotron loss time t{p) = 'dm'lc^ / {ArlB\p) , where me 
is the electron mass, ro « 2.8 x 10~^'^cm is the classi- 
cal electron radius, c denotes the speed of light, and Bd 
is the effective magnetic field in the downstream region. 
For given frequency v the momentum p of the radiating 
electron is approximately given as p cx v / B^. 

In the limit of a strongly fluctuating magnetic field 
around the shock on all scales, we also assume that the 
diffusion coefficient k{p) is given by the Bohm limit, 
^{v) = Pbi^/S, where pB and v are the particle gyrora- 
dius and velocity, respectively. Therefore I2 is determined 
by i?d and v which can be conveniently turned into an 
equation for in terms of I2 and v (Berezhko & Volk 
2004a): 



Pd = [^myi{Aerlll)Y/\VY^ -5)-^'\ 
with 



(2) 



(3) 



where e is the proton charge. These are the main relations 

which wc arc going to use in order to determine the am- 
plified field. The Chandra observations yield J,y{e^,p) in 
terms of the angular distance Atp = {p — Rs)/d which can 
be converted to p for known source distance d. Given the 
Ji,(/9) -profiles for a given X-ray energy we shall then 
fit the data to the model curve Eq.(l) to determine I2, 
and thus Pd- In some of the published observations, men- 
tioned in the Introduction, only the exponential widths of 
the profiles are quoted. In those cases we shall approxi- 
mately equate L with this width. 

2.1. Tycho's SNR 

In the case of Tycho's SNR the circumstellar gas density is 
unexpectedly inhomogeneous for a type la SN. This results 
in significant azimuthal deviations from spherical symme- 
try both in shape and expansion velocity of the SNR shock 
(e.g. Reynoso et al. 1999), even though to first approxi- 
mation the remnant can be considered as spherical. There 
are also indications that the local radio spectral index a, 
derived from the ratio of the emissions at A = 20 cm and 
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A = 90, is nonuniform (Katz-Stone et al. 2000) and in 

several locations quite low, near a = 0.5. The well-defined 
spatially integrated radio spectrum on the other hand is 
rather steep (a w 0.61) if fitted by a power law, and indi- 
cates curvature (Reynolds & Ellison 1992). It is therefore 
perhaps not surprising that also the radial brightness pro- 
files across the shock are not identical in different regions 
of the shock circumference. Two such profiles were pre- 
sented by Hwang et al. (2002). 

There are additional reasons for different widths of 
such profiles apart from inhomogeneities of the circumstel- 
lar density. Being an example of a type la SN, the low-mass 
progenitor star should not have had a significant mass loss 
in the form of a wind with a Parker spiral-type circum- 
stellar field. To the extent that the density of the environ- 
ment is uniform, also the ambient magnetic field should 
be uniform, ideally resulting in a dipolar topology of the 
distribution of accelerated nuclear particles (Volk 1997; 
Volk et al. 2003). This should also lead to a correspond- 
ing distribution of field amplification in azimuth and thus 
to different spatial scales of the radial synchrotron pro- 
files. Even though such a simple geometry is not realized 
in Tycho's SNR, the underlying physical processes should 
nevertheless lead to azimuthal variations of the radial syn- 
chrotron emissivity scales in a quasi-statistical way. 

As indicated in the Introduction a global value of the 
effective interior magnetic field - and thus the sharpness of 
the radial profile at the shock - can be independently de- 
termined from a comparison of the volume-integrated the- 
oretical synchrotron spectrum with the one observed. In 
fact, the nonlinearity of the diffusive acceleration process 
for CR nucleons leads to a spectral steepening and even 
to a curvature of the energy spectrum of accelerated elec- 
trons for rigidities corresponding to sub/trans-relativistic 
protons. This well-known feature (e.g. Ellison & Eichler 
1984; Drury 1983; Volk 1984; Blandford & Eichler 1987; 
Reynolds & Ellison 1992; Berezhko et al. 1996) implies 
a steeper and ciirved synchrotron spectrum in the radio 
range compared to the test particle limit, if the radiating 
electrons are indeed in the above low rigidity range. An 
observed integrated spectrum with this feature then im- 
plies a minimum magnetic field strength. But this is only 
the most obvious nonlinear feature. In fact, the entire fre- 
quency range of the synchrotron spectrum must be con- 
sistent with this magnetic field strength, especially also in 
the hard X-ray region. This fixes both the amount of shock 
modification, i.e. the proton injection rate in the theory, 
and the global effective magnetic field strength (Berezhko 
et al. 2002). It is remarkable that with these three param- 
eters it is possible to describe the measurements of the 
entire synchrotron spectrum, which in the case of Cas A 
include also the mm and mid-infrared range. 

As a consequence we have two experimentally indepen- 
dent methods to determine the effective magnetic field, 
a local one and a global one. Several independent radial 
profiles even yield several independent local values as a 
function of azimuth. We shall in the following investigate 
the consistency of these values for Tycho's SNR. 



The radial width I2 of the three-dimensional emissivity 
distribution is smaller by a factor of about 7 than the 
radial brightness scale of the optically thin emission. It 
is simple to verify that for the considered X-ray energies, 
which lie in the cutoff region of the integrated synchrotron 
brightness, the case of strong losses (Volk et al. 1981) 
is relevant, where the emissivity width is given by the 
Umiting case of Eqs.(2, 3) (Berezhko & Volk 2004a) 

h = \/i^iTi, (4) 

i.e. by the postshock diffusion length in a synchrotron loss 
time. In the Bohm diffusion limit for the accelerating par- 
ticles in Tycho's SNR (which at the present epoch is in 
the transition between the free expansion phase and the 
Sedov phase), I2 is then given by 

h = [3my/{4erlBl)]'/'. (5) 

In other words, the radial emissivity profile is 

parametrized by alone in this limit. 

We shall start with the Chandra data that were given 
in Fig. 4 of Hwang et al. (2002) in the form of radial pro- 
files of the 4-6 keV continuum brightness for two azimuthal 
sectors in the northwest (NW) and the southwest (SW). 
This energy range excludes most of the emission lines. The 
profiles show a pronounced outer rim, presumably at the 
location of the remnant's forward shock. 

In the lowest approximation the narrower rim in the 
NW has an exponential width of about 3.5 arcsec, whereas 
the SW rim extends over at least 4.7 arcsec. This cor- 
responds to about 1.6 and 2.1 percent of projected ra- 
dius. Adopting a distance of d = 2.3 kpc and using 
I2 ^ L/7, this implies I2 ^ (1.7 ± 0.7) x 10^'^ cm and 
I2 w (2.3 ± 0.7) X 10^6 cm for the NW and SW profile, 
respectively. Using the full Eqs. (2) and (3) with Vg = 
3100 km/s, and a shock compression ratio cr = 6, from 
Volk et al. (2002), and taking = hu ^ 4: keV, we obtain 
Ed ~ 404 (+169 - 82) and ~ 332 (+88 - 53) /iG, 
for the two respective sectors. The shock sectors are obvi- 
ously independent except possibly the innermost regions 
near the center which play no role in the present argument. 
Nevertheless one can conclude that the magnetic field val- 
ues within these regions are roughly the same within the 
errors. 

Bamba (2004) has given an exponential width L = 
0.073(-|-0.010- 0.009) pc for a profile in the NW, in the 2 
- 10 keV band. Using e,y = 2 keV and I2 = L/7 we obtain 
from Eqs. (2) and (3) B^ = 301(+28 - 25) /xG. 

We have also fitted both the widest and most narrow 
template distributions, cf. Eqs. (1-3), to the Hwang et al. 
profiles. To do this we assumed a Gaussian distribution 
of the statistical errors in the measurement of the bright- 
ness distribution , i.e. the errors to be given by the 
square root of the discrete values given in the paper by 
Hwang et al. (2002), together with a positional accuracy 
of the ACIS array of Chandra of 0.5 arcsec. The results 
for the two radial profiles are given in Fig. 1. This im- 
plies Bd = 331(+182 - 94) /xG for the SW profile, and 
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Fig. 1. The Chandra 4 — 6 keV continuum brightness radial 

profile of Tycho's SNR observed by Hwang et al. (2002) . Upper 
panel (a): southwest azimuthal sector. Lower panel (b): north- 
west azimuthal sector. The statistical error in the brightness 
values was assumed to equal the square root of the values, 
whereas for the angular position an error of 0.5" was assumed. 
The data are fitted to the projection of the exponential emissiv- 
ity profile described by the model of Eq. (1) in the downstream 
region (solid line). The fit has a = 17.3, with /dof — 1.73 
and = 6.56, with x? /dof = 0.82 for (a) and (b) respectively. 
The dotted lines indicate the la deviations. 

Ba = 272(+117 - 70) iiG for the profile in the NW. Note 

that the vahics jdo) = 1.73 and /dof = 0.82 for the 
cases represented in Fig. la and lb respectively show that 
the data are fitted by the formula (1) quite well. Again 
one can see that within the uncertainties the field values 
in these two region are the same. 

For the global determination of the effective down- 
stream field Bi\ wc compare the theoretical synchrotron 
spectrum with the observed spatially integrated spectrum 
(see Fig. 2). To explore the possible bandwidth of Ba- 
values we approximate the data with a slightly softer ra- 
dio spectrum than adopted by Volk et al. (2002). In the 
framework of the strongly nonlinear system of equations, 
the required stronger shock modification by the nuclear 
component is achieved by a somewhat increased proton in- 
jection rate. The theoretically implied increase of the mag- 
netic field from 240 to 360 then requires a reduction 
of the amplitude of the electron momentum distribution 
since the radio electrons are not subject to synchrotron 




I I I I I I I I I \ I \ d 

8 10 12 14 16 18 

og f, Hz 

Fig. 2. Synchrotron spectral energy density vSi, of Tycho's 
SNR as a function of frequency v. The solid curve corresponds 

to a fit with a downstream magnetic field of = 240 /iG (Volk 
et al. 2002) which slightly overestimates the X-ray continuum 
at 1/ > 10^* Hz and therefore underestimates the field. The 
more strongly nonlinearly modified case, with slightly steeper 
radio spectrum ( dashed curve) requires a higher B-field of Ba = 
360 /iG but somewhat underestimates the X-ray brightness. 
The actual B-field value lies in between these cases. 

cooling at the present epoch. This leads to a decreased 
electron/proton ratio in the energy range where radia- 
tive cooling is unimportant. Most importantly however, 
the cooling region of the synchrotron emission spectrum 
V > 10^^ Hz will be lowered in amplitude by this decreas- 
ing electron/proton ratio. 

Taking into account, that the electron distribution 
function is roughly a power law /e a we can approxi- 
mately write 

cx p'fMBdh (6) 

for V PS {p/mc)'^eBd/mc, where fc{p) is the spatially aver- 
aged electron distribution in the remnant interior. The low 
frequency part of the spectrum S^{p < 10^^ Hz) is pro- 
duced by electrons which do not suffer synchrotron losses. 
Therefore I2 ~ 0.1i?s and Si, cx r]Kcpiy~"B2'^^, where rj 
is proton injection rate, Kep is the electron-to-proton ra- 
tio, and a = {q — 3)/2. To have an equally good fit for 
the radio data the electron-to-proton ratio Kcp should be 
reduced for a higher downstream field value B^- 

At higher frequencies u > 10^^ Hz the synchrotron 
emission is produced by electrons which suffer signif- 
icant synchrotron losses. For v < fmax = 10^^ Hz 
I2 ~ U2T2, which corresponds to the so-called weak loss 
case (Berezhko & Volk 2004a). Since T2 oc jr^B^^ we 
have Si, oc r]KepP^~'^ / B^ oc r]KepV~°'B'^~^ , where a = 
{q — 2)/2. Since the power law index q is close to 4, 
Si, oc r]Kepy~^, independently of B^- The reduction of 
the electron- to-proton ratio ii'cp oc iJ^''"^^/^ for higher 
field i?d leads to a decrease of the overall synchrotron flux 
5, (see Fig. 2). 
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Fig. 3. Comparison of the low-field Bd = 240 /iG (solid curve) 
and high- field -Bd = 360 nG [dashed curve) radial profiles of 
the full numerical solution in spherical symmetry with the ob- 
servational Chandra profiles (a) and (b) from Fig. 1. 



The reduction of electron production, discussed above, 
leads to a proportional reduction of the synchrotron emis- 
sion in the cooling range as shown by the dashed curve 
in Fig. 2 and implies an underestimate of the X-ray syn- 
chrotron emission compared to the slight overestimate in 
Volk et al. (2002, see Figs. 3 and 5). The spectral fit in 
that earlier paper had on purpose not been fine-tuned, be- 
cause the aim had been to demonstrate field ampHfication 
in principle. The optimum magnetic field strength then 
lies between the two cases, roughly at 300 ^G. This leads 
to an optimum value w 300 ± 60 /uG from acceleration 
theory. 

The two extremes of the globally determined magnetic 

field from Fig. 2 are then compared to the local Chandra 
data in Fig. 3 by plotting the full numerical solutions for 
the SNR morphology, as it results from the combination 
of gas dynamics and time-dependent kinetic theory, to- 
gether with the data from Fig. 1. The agreement of global 
morphology and local profiles is reasonable even though 
the differences between the two profiles show that the de- 
viations from spherical symmetry are significant. Such a 
comparison is of course incomplete since only two radial 
profiles could be used. Further Chandra data are needed to 
evaluate in detail the degree of agreement between theory 
and experiment. 



Fig. 4. Chandra ACIS spectroscopic array image of the north- 
ern and northwestern part of Tycho's SNR in the 4-6 keV 
band. Shown are 6 selected boundary regions (boxes) that are 
assumed to contain portions of the shock. Each region defines 
a radial profile from which a value of Bd is derived. The sum 
of all 6 profiles defines an average profile which is shown in 
Fig. 5. 



Finally, we have investigated Chandra X-ray images 
in the north and northwest, obtained for X-ray energies 
4 < e,y < 6 keV with the ACIS spectroscopic array 
from the Chandra archive, using the morphology tem- 
plate. This corresponds to 6 profiles through the remnant's 
outer boimdary as indicated in the imsmoothcd broad- 
band Chandra ACIS image in Fig. 4, showing a partial 
image of the SNR. In addition we show the co-added data 
from the 6 profiles. The corresponding data are fitted to 
the template. The individual results for vary between 
150 /xG and 373 /xG, with typical erros of (- 34 % + 70 
%). Since the number of events per bin is partly small, 
of the order of a few in the outer fianks of the individ- 
ual distributions, we have also coadded the profiles to ob- 
tain a sum profile. It is shown in Fig. 5. Even though the 
individual low-statistics profiles lead to a sizable scatter 
for Bd, the average profile determines a mean value of 
Bd = 273(+49 - 37) /iC with a very high fitting qual- 
ity which agrees rather well with the optimized field value 
from acceleration theory. 

Note, that the actual shock position relative to the ob- 
served brightness profile is not known. Therefore it is used 
in the fitting procedure as a free parameter. As a result of 
the best fit the shock position corresponds to Atp = 0.3 
and 0.6 arcsec for the cases presented in Fig. la and lb 
respectively, and Aip = — 1.9 arcsec for the case, presented 
in Fig. 5. The overall situation is then as follows: the sim- 
ple approximation that determines an exponential width 
L and then uses Eqs.(2) and (3) to determine I2, yields B^- 
values that agree well with each other and with the global 
spectral determination. The same is true for the fit of the 
profiles from the archive to the template given by Eq.(l). 
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Fig. 5. Sum profile of the 6 radial X-ray profiles imaged in 
Fig. 4. The data are fitted to the projection of the exponen- 
tial emissivity profile described by the model of Eq. 3 in the 

downstream region {solid line). The fit has a = 9.73, with 
/dof = 0.69. The dotted lines indicate the la deviations. 



Therefore we can also conclude that the coincidence of 
the postshock and the global field values confirni that the 
field amplification during the SNR evolution scales like 

B"^ = canst x . 

2.1.1. Gamma-ray prediction 

For completeness we shall also give the resulting prediction 
for the total integral gamma-ray energy flux from Tycho's 
SNR for the above extreme values of B^, equal to 240 /xG 
and 360 /xG, respectively, at the same gas density (Fig. 6). 
Compared to the earlier prediction (Volk et al. 2002), the 
hadronic 7r°-decay flux changes very little, as expected, 
whereas the Inverse Compton flux changes by a factor of 
about 7 between 240 yuG and 360 ^G in effective magnetic 
field strength, given the observed synchrotron flux. Thus a 
detection of Tycho's SNR in TeV gamma rays close to the 
predicted hadronic flux would again indicate a hadronic 
CR source - as the entire discussion in terms of acceler- 
ation theory implies in any case. However, the hadronic 
gamma-ray emission is quite sensitive to the ambicint gas 
density which may not have been fixed sufficiently well 
yet from radio observations (e.g. Rcynoso et al. 1997). We 
note however that the nominal value for the mean ambient 
hydrogen density iVn = 0.5 cm~^, chosen in Fig. 6, is on 
the low rather than on the high side. 

2.2. RCW 86 

This SNR, more generically designated as G315.4-2.3, is 

a shell-type SNR of large angular size ^ 45 arcmin, with 
a nonthermal radio spectrum, at an apparent kinematic 
distance oi d = 2.8 kpc found from optical observa- 
tions (Rosado et al. 1996). The distance might in fact be 
smaller, and might be as low as d ~ 1 kpc (see e.g. Dickel 
et al. 2001; Borkowskiet al. 2001, for more recent reviews). 



HA 



W \ 

^H-CT 




Fig. 6. Predicted 7-ray spectral energy distribution, as a func- 
tion of 7-ray energy: IC emission {dashed lines), Nonthermal 
Bremsstrahlung (NB, dash-dotted lines), and 7r"-decay {solid 
lines). The observed 3o- 7-ray upper limits (W - Whipple 
(Buckley et al. 1998), H-CT - HEGRA lACT system 
(Aharonian et al. 2001b)), and the 95% confidence HA - 
HEGRA AIROBICC upper limit (Prahl et al. 1997) are shown 
together with the detection/upper limit claimed by the Shalon 
group (Sinitsyna et al. 2003). Thick lines correspond to a down- 
stream field of -Bd = 240 /itG, whereas thin lines stem from 
assuming the maximum field value Bd = 360 /xG. 



The designation RCW 86 is an optical one and refers to a 
bright complex of optical emission filaments in the south- 
west part of the more or less circular and complete shell. 
This southwest part also exhibits the brightest radio and 
X-ray emission. G315.4-2.3 has often been compared to 
Tycho's SNR. However, it lacks the sharp outer rims that 
delineate much of Tycho's SNR silhouette (cf. Dickel et 
al. 2001), in particular also in hard X-rays (Fig. 4). 

The X-ray observations above 2 keV are well described 
by a dominant synchrotron continuum and a strong iron 
Ka line. Using Chandra, Rho et al. (2002) in particular 
investigated spatial profiles of the hard (2 < Cj, < 8 keV) 
X-ray emission in the RCW 86 region proper, one of which 
we consider appropriate to represent a strong shock since 
it has all the expected properties: a sharp rise in hard X- 
rays - accompanied by a much smoother rise in the 1.4 
GHz radio emission - followed by a somewhat slower de- 
cline towards the same level of background. For this region 
in the west of RCW 86 the 1.4 MHz emission profile has a 
more than ten times larger decay scale than the X-ray pro- 
file and may well be strongly influenced by adiabatic ex- 
pansion effects, masking the comparatively weak radiative 
losses at those electron energies. We shall concentrate on 
this single profile which in any case yields the largest value 
of the effective field available for this source. The data for 
this profile appear to also have been used by Bamba et al. 
(2003b) and Bamba (2004) in their analysis. Their expo- 
nential angular width leads to L = 0.29 pc for the sharpest 
profile, that is very close to our estimate L = 0.26 pc 
from the Rho et al. (2002) profile. However, even this pro- 
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Fig. 7. Fit of the data (histogram) corresponding to the Rho 
ct al. (2002) profile from Chandra in the 2-8 keV X-ray band 
for RCW 86 to the template cf. Eq.(l) (solid line). The la 
deviations are given by the dotted curves. 



file is in the interior of RCW 86, not at the rim. It is 
possible that this part of the shock surface has a cur- 
vature due to an inhomogcnuity of surrounding medium 
which is quite different from that of an ideal spherical sur- 
face. We therefore have to consider the likelihood of non- 
sphcrical geometry and thus of incomplete dc-projection 
with a ratio L/I2 < 7, as in the case of RX J1713. 7-3946 
below. Using nevertheless Eq.(2) for the time being, we 
obtain « 108 /iG for our estimate from the Rho ct 
al. data and 100 /xG for the Bamba et al. width. 

The shock has been assumed to move with a velocity of 
Vs = 800 km/sec into a low-density medium with hydro- 
gen density A^h = 0.3 cm~'^ at a kinematic distance of 
d = S kpc, cf. Borkowski et al. (2001). 

If we apply the model fit cf. Eq.(l) to the top profile of 
Fig. 8 of Rho et al. (2002) we obtain B^ « 99(+46-26) fiG 
with a x^/dof = 0.35, using 25 dof and assuming spherical 
geometry as our preferred value for full de-projection (see 
Fig. 7). 

If we were to assume a distance of d = I kpc to 
RCW 86 instead, the linear scales would decrease by a fac- 
tor of about 3, increasing the effective field by a factor of 
about 2 to iJd = 231 ^G. Whereas a value of B^ = 100 /xG 
corresponds to a comparatively modest field amplification 
of about a factor of 5, we must therefore count on a very 
strong amplification in the case that the source is nearby. 
This would be no surprise, assuming that for the close dis- 
tance the source might be as young as 2000 yr, instead of 
~ 10^ yr for the large distance (Borkowski et al. 2001), 
resulting in a high shock velocity and correspondingly ex- 
pected Pc ^ pV^. 

In reality however we must count on somewhat larger 
scales I2 of the emissivity for the observed brightness scale 
L. A ratio I2/L = 1/4 or even I2/L = 1/2 would reduce the 
above field amplification by a factor 0.69 or even 0.43, re- 
spectively. These appear realistic correction factors, given 



the fact that the shock transition lies in the outer part of 
the shell. 



2.3. Kepler's SNR 

Kepler's SNR is 400 yr old and has an angular radius of « 
100 arcsec at a distance d ~ 4.8±1.4 kpc (Reynoso & Goss 
1999). This implies a linear radius Rg ~ 2.3±0.7 pc, some- 
what smaller than Tycho's SNR. Like the optical and the 
radio continuum emission the X-ray brightness distribu- 
tion from Kepler's SNR is rather asymmetric with stronger 
emission in the north and northwest (e.g. DcLancy ct 
al. 2002; Cassam-Chenai' et al. 2004a). The most obvi- 
ous nonthermal hard X-ray continuum emission in the 4 - 
6 keV range, as observed with XMM-Newton, is found in 
a southeastern region, rather at the opposite side of the 
remnant (Cassam-Chenaz et al. 2004a), where also thin fil- 
aments have been observed with Chandra. The spectra of 
these filaments have been found to be hard without line- 
structures by Bamba et al. (2003b) and by Bamba (2004), 
from Chandra archival data, with the conclusion that the 
filaments are the result of synchrotron emission. One pro- 
file has again the characteristics of a shock transition, for 
which they show a scale length of about L = 0.11 pc in 
the 4.0-10.0 keV band. 

Applying our de-projection model to this source de- 
spite its asymmetry, we obtain I2 ~ 1.6 x 10~^ pc. Cassam- 
Chenaz et al. (2004a) identify this southeastern part of the 
SNR as a very low density medium with A^h < 0.15 cm~^ 
and a shock velocity Vs ~ 4550 km/s. Using this velocity 
in Eq.(3) for 6, and assuming an overall shock compres- 
sion ratio of CT = 6, we obtain from Eq.(2) the substantial 
field Bd « 215 /iG, with an estimated error of less than 
20%. Thus we conclude that there is substantial field am- 
plification in Kepler's SNR. 



2.4. SNR RX J1713.7-3946 

This very large « 1 degree diamciter shell-type SNR in 
the Galactic plane was discovered in the ROSAT all-sky 
survey (Pfeffermann & Aschenbach 1996). The hard X-ray 
emission detected with ASCA turned out to be purely non- 
thermal and did not reveal any thermal emission (Koyama 
et al. 1997; Slane et al. 1999). Recently a weak thermal 
component may have been detected in the interior of the 
remnant using large field-of-view X-ray instruments, in- 
cluding the RXTE Proportional Counter Array (PGA) 
(Pannuti et al. 2003). The CANGAROO experiment re- 
ported a detection in TeV 7-rays at a level of 70% of the 
Crab Nebula (Muraishi et al. 2000; Enomoto et al. 2002). 

The H.E.S.S. experiment has confirmed (Aharonian et 
al. 2004) the CANGAROO detection. The gamma-ray and 
X-ray images have similar shapes. There appears to be no 
doubt about the reality of this source in TeV gamma-rays. 

The SNR is complex. New CO-observations by Fukui 
et al. (2003) suggest the source to be rather close, at a 
distance of about d = 1 kpc. A short distance of d = 
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1.3 ± 0.4 kpc has recently also been inferred by Cassam- 

Chenai' (2004b) from an a detailed investigation of the 
X-ray absorbing material in the remnant together with 
the CO-observations. 

Hard X-ray (1 < f,y < 5 kcV) observations with 
Chandra by Uchiyama et al. (2003) and Lazendic et al. 
(2004) revealed filaments in the remnant interior. They 
appear embedded in a diffuse plateau emission (Cassam- 
Chenai' 2004b). According to Uchiyama et al. the most 
prominent and largest filament has a very small width of 
about 20 arcsec. For a distance of = 1 kpc this translates 
into a spatial scale of L = 0.1 pc. At a different portion 
of this filament Lazendic et al. (2004) inferred a width of 
about 40 arcsec, i.e. twice as large. We note however that 
the width drived from the observed profiles depends ap- 
preciably on the assumed background level. For a lower 
background level than that assumed by Uchiyama et al. 
(2003) and Lazendic et al. (2004) the experimental width 
may be considerably higher, with a corresponding mag- 
netic field value that is considerably lower than estimated 
below. 

It is likely but not certain that the filament measured 
corresponds to a shock. A tangential discontinuity due to 
an internal structure of the circumstcllar medium before 
the explosion or a discontinuity separating shocked ejecta 
from the shocked circumstellar medium cannot be ex- 
cluded. To this physical possibility one must add the prob- 
ability of shear- type three-dimensional convective flows, 
distinct from spherical expansion. In particular, the pro- 
jection of such a structure along the line of sight might 
be far different in scale from that for a spherical structure 
with a radius comparable to the dimension of the overall 
object. It is indeed quite possible that the Chandra pro- 
file cuts have not crossed the outer shock, and therefore 
do not necessarily delineate the main field amplification 
region caused by the accelerated particles. 

Nevertheless, we can attempt some estimates for dif- 
ferent possible shock situations regarding the observed fil- 
ament: (i) a shock transition of scale L = 0.1 pc, assuming 
a compression ratio cr = 5 and a speed of = 4000 km/s, 
and a low upstream density A^h = 2 x 10~^ cm~^ (Cassam- 
Chenai' 2004b) for an X-ray energy e,^ = 1 kcV gives a 
maximum « 271 /xG; it reduces to iJa « 74 jiG with- 
out any de-projection, I2 = L, and to ^ 118 /xG for 
I2 = L/2; (ii) for a shock in a denser medium A'h ~ 
1.3 cm~"^, with velocity Vg = 1000 km/s, the maximum 
Bd = 213 fxG; it reduces to B^ 58 /xG for I2 = L, and to 
i?d = 94 /iG for I2 = L/2. For the twice thicker Lazendic 
et al. filament the foregoing field values would reduce by 
a factor of about 0.6. This range of field strengths might 
provide an idea what to expect from interior structures in- 
side this SNR. New high-resolution measurements appear 
necessary for the global picture. 



2.5. Cas A 

Extending Bamba et al. (2003b), Bamba (2004) gives two 
Chandra profiles through the southeastern rim at X-ray 
energies 5 < < 10 keV whose average width is L w 
0.028 pc and which we convert to Z2 ~ 4 x 10~^ pc. Using 
the hydrodynamic solution of Berezhko et al. (2003b) for 
the remnant dynamics, with a shock compression ratio 
cr « 6 and a shock speed Vg = 2400 km/s, we obtain 
Bd ~ 485 /xG, taking = 5 kcV. If we were to take 
Vs = 3000 km/s we would obtain B^ ~ 496 /iG. 

Vink & Laming (2003) have analyzed a profile with 
L « 0.025 pc through the northeastern rim at 4 < < 
6 kcV, and assumed a shock velocity of = 5000 km/s. 
For Vs = 3000 km/s and using this scale L, we calculate 
Ba w 550 fiG (Berezhko & Volk 2004a). 

These values from the two different profiles agree quite 
well with each other. They differ by less than 15 % 
from the spectrally determined value B^ ~ 480 /xG from 
Berezhko et al. (2003b). 

2.6. SN 1006 

For SN 1006 Bamba et al. (2003a, 2004) for > 2 keV 
give exponential widths for the sharpest rim profile L = 
19 arcsec, which for the distance d = 2.2 kpc translates 
into L = 0.2 pc. Using a shock velocity Vs = 3200 km/s 
and a shock compression ratio a = 6.3, Eq.(2) yields a 
downstream field value B^ « 143 /xG. 

According to Bamba et al. (2003a, 2004) and Bamba 
(2004) there is an even sharper profile with L = 0.126 pc, 
which is however not a rim profile. With I2 = L/7 this 
leads to « 195 /uG. 

These numbers can be compared to the average field 
value 120 /xG predicted from the overall synchrotron spec- 
trum by Berezhko et al. (2002), and depending very 
weakly on the ambient gas density. In fact, reconsider- 
ing the spectral determination of B^, wo can achieve a 
preferred field value of 160 /xG (Ksenofontov et al. 2004). 

Taking into account the data uncertainty, which is 
about 25%, we can conclude that the above numbers de- 
rived from the Bamba et al. (2003a, 2004) data are con- 
sistent with the values derived from the synchrotron spec- 
trum. Together they confirm the overall field amplification 
picture for SN 1006. On the other hand, there may also 
be a physical reason for the differences in the values of B^ 
across the shock surface. This might occur if in the po- 
lar caps independent magnetic flux tubes exist, between 
which the field strength and probably also the generating 
CR pressure Pc can fluctuate considerably. 

3. Discussion 

All four newly calculated sources show amplified fields. 

In addition Tycho's SNR furnishes a third case where 
the field strength B^ determined from small-scale syn- 
chrotron features agrees with that predicted by comparing 
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the global synchrotron spectrum with nonlinear accelera- 
tion theory. 

Although the field amplifications may be relatively 
small for RCW 86 and RX J1713. 7-3946 and arc difficult 
to pin down, all six young Galactic SNRs with known fil- 
amentary scales exhibit field amplification. There appears 
to be no exception to this "rule" up to now. 

As an important prerequisite the characteristic nonlin- 
ear emission features require a dominant nuclear compo- 
nent of accelerated particles. At the moment this conclu- 
sion is in part still a theoretical one: although the 7-ray 
fiux at 1 TeV from Gas A, detected by the HEGRA ex- 
periment (Aharonian et al. 2001a), agrees with the theo- 
retically calculated hadronic flux there is as yet no con- 
firmation by another experiment. SN 1006, reported by 
the CANGAROO experiment (Tanimori et al. 1998, 2001) 
could not be detected by the H.E.S.S. experiment as a 
TeV source until now (Masterson et al. 2003; Rowell et 
al. 2004). An inverse Compton 7-ray emission scenario 
in a low magnetic field of order w 10 /iG (e.g. Pohl 
1996, Mastichiadis & de Jager 1996; Tanimori et al. 1998, 
2001) is clearly excluded by the H.E.S.S. upper limits. 
The H.E.S.S. upper limit is almost one order of magni- 
tude lower than the published CANGAROO fiux. Given 
the range of ambient densities 0.05 < A^h/(1 cm~^) < 0.3 
it suggests a low density ISM with TVh < 0.1 cm~^ 
(Ksenofontov et al. 2004) instead of the high value Nu = 
0.3 cm"^, which happened to fit the CANGAROO data 
(Berezhko et al. 2002). Finally, the predicted hadronic 
gamma-ray fiux for Tycho's SNR (Volk et al. 2002) is still 
below the upper limit measured by HEGRA (Aharonian 
et al. 2001b). Only RX J1713. 7-3946 has been detected by 
two independent gamma-ray experiments. 

In all cases the gamma-ray emission depends strongly 
on "external" astronomical parameters, like the ambi- 
ent gas density, and therefore requires extensive multi- 
wavelength investigations with significant inherent uncer- 
tainties. 

Regarding the magnetic field strengths derived we can 
attempt to learn something from these results by empiri- 
cally correlating the magnetic field pressure B^/ (87r) with 
the ram pressure poV^ of the upstream gas in the shock 
frame of reference. Assuming that Bj/(87r) reaches some 
fraction of the CR pressure Pc at the shock, and assum- 
ing very efficient acceleration so that Pc = epoV^, with 
e = 0(1), we expect a constant ratio B'^/ {%-k pqV^) as 
a function of shock velocity Vs for an individual object. 
This appears to us to be the physically most natural re- 
sult, given the overall dynamics in the shock (see also the 
discussion in Volk et al. 2002). Berezhko & Volk (2004b) 
have used this assumption to calculate the proton and 
electron spectra for a generic SNR with field amplifica- 
tion. The result is that the proton component reaches the 
so-called knee in the Galactic CR spectrum with a maxi- 
mum momentum p 

max 

10*^m„c. 



In numbers the ratio of magnetic and kinetic energy 
densities is given by: 



1.7x10" 



[I?a/(100 ijG)f 



[Nnl{l cm-3)][v;/(103 km/s) 



,(7) 



where B^ may itself implicitly depend on Vg if it is deter- 
mined from Eqs.(2) and (3). 

On the other hand. Bell (2004) has argued for a pro- 
portionality B\/{Si-k) oc {ys/c)pQV^ from his quasi-MHD 
simulations of the nonrcsonant right-hand polarized low 
frequency MHD mode. This implies a linear increase of 
-B^/(87rpoK^) with V^. Bell has suggested that accelera- 
tion even considerably beyond the knec> may be possible in 
this way in very young SNRs expanding into dense circum- 
stellar material. This extends earlier arguments put for- 
ward by Bell & Lucek (2001) who estimated the maximum 
achievable energy cpmax- However, the results of Berezhko 
& Volk (2004b) indicate that the energy content of the 
fastest part of the supernova ejecta, which produces CRs 
with energies above the knee energy 3 x 10^'"' eV, is so small 
that the resultant CR spectrum produced in SNRs is very 
steep N oc e~^. Therefore these highest energy CRs play 
no signigifant role in the formation of the Galactic CR 
spectrum above the knee. The reason is the shortness of 
the time period during the initial SNR phase over which 
a very high energy fiux density pV^/2 into the shock is 
maintained. 

Obviously we can not follow an individual source as it 
evolves in time. The only possibility is to make the strong 
assumption that different objects at different epochs of 
their expansion history behave like a single source at such 
points in its evolution. I.e. we can substitute an ensem- 
ble average for a time average. Given that we consider 
here even different types of SN explosions (Type la for 
SN 1006, Tycho, and possibly Kepler, probably type II 
for RX J1713.7-3946 and RCW 86, and probably type lb 
for Gas A), conclusions must be drawn with this fact in 
mind. 

Fig. 8 shows a plot of B'^/{8iTpoV^) vs. Vs for the 
six sources discussed in the previous section. For Tycho's 
SNR, SN 1006, and Cas A the available globally deter- 
mined field values were used. For the sources with rim 
profiles (Tycho's SNR, SN 1006, Kepler's SNR, and Cas 
A) we have assumed the full de-projection L = 7/2, and 
ATh = 0.1 cm-3 and A^h = 2 cm'^ for SN 1006 and Cas A, 
respectively. 

For RCW 86 we have assumed a range of uncertainty 
7/4 < 7I2/L < 7/2 (the four points in Fig. 8 correspond 
to L/I2 = 1,2,4 and 7, respectively), and for RX J1713.7- 
3946 we have taken the even more pessimistic range 7/2 < 
7I2/L < 7 (the three points in Fig. 8 correspond to L/I2 = 
1, 2 and 7). These values correspond to the discussion of 
these two sources in the previous section. 

For RX J1713.7-3946 we have in addition plotted the 
two alternatives corresponding to = 1000 km/s and 
A^H « 1.3 cm~3 on the one hand, and 14 = 4000 km/s 
and Nu w 2 x 10~^ cm"^ on the other. 
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Fig. 8. Ratio of downstream magnetic field pressure _Bj/(87r) 
to preshock gas ram pressure poV^ (in percent) vs. shock veloc- 
ity (in units of 10^ km/s) for the six young SNRs that show 
field amphfication. For the objects RCW 86 and RX J1713.7- 
3946 a range of normalized magnetic field pressures is given. 
(For RCW 86 the four points correspond to L/h = 1, 2, 4 and 
7, respectively, and the solid vertical line segment corresponds 
to the assumed range of uncertainty 7/4 < 7I2/L < 7/2. For 
RX J1713. 7-3946 the three points correspond to L/h = 1, 2 
and 7, respectively, and the solid vertical line segment corre- 
sponds to the uncertainty range 7/2 < Yh/L < 7) For this 
remnant two alternative shock states are depicted: the low ve- 
locity case (a) corresponds rather to the cloudy medium in 
the north and west of the remnant. The high velocity case (b) 
is more representative for the very low density regions in the 
southern projection of the remnant, the fully projected case 
L/I2 = 7 corresponding to a value of 40, out of scale in the 
figure. 



The present uncertainties in the data, which are not 
more accurate than 25% and cover only a rather limited 
range of the shock speed Vg, do not allow a distinction 
between a constant ratio B^/{8TrpoV^) and one with a 
linear growth of B'^/{87rpoV^) with Vs. However, the data 
are highly likely to exclude a decrease of B^/{8TrpoVg ) 
with increasing Vg. Obviously this means that with grow- 
ing age, that is decreasing shock velocity, the field am- 
plification in SNRs decreases about linearly with the 
shock velocity. Fig. 8 indicates that in all young SNRs 
with large enough shock speed Vs > 10'^ km/s the ef- 
fective (amplified) magnetic field energy density is near 
Bl/iSn) = 3.5 X IQ-^poV^^, the value for Kepler's SNR 
probably having the largest uncertainty. It is important 
to note that in all four of the best defined cases (SN 1006, 
Gas A, Tycho's SNR and Kepler's SNR), which are char- 



acterized by rather close values of the shock speed, the 
normalized magnetic field pressures, extracted from the 
data, are also close to each other. The large uncertain- 
ties in the determined magnetic field values for the cases 
of RCW 86 and RX J1713.7-3946 are due to much more 
poorly known values of their relevant parameters, such as 
age, distance, shock speed and ISM gas density. 

We thus tentatively conclude that field amplification 
is confined to young remnants. This also implies that the 
escape of the highest energy nuclear particles, accelerated 
at an earlier epoch with high effective field, becomes pro- 
gressively important as the remnant age increases. This 
implies a lower and lower cutoff of the resulting gamma- 
ray spectrum from a SNR with age. As a result, the high- 
est energy nuclear particles are to be found in SNRs that 
have just reached the Sedov phase. Older remnants are 
increasingly unable to confine them. 

The amplified magnetic field also leads to a depression 
of the Inverse Compton and nonthermal Bremsstrahlung 
gamma-ray emission relative to the hadronic emission 
for young remnants, given the synchrotron emission. 
Nevertheless the hadronic gamma-ray emission depends 
strongly on iVn, essentially cx (Volk 2003). Thus a 
strong source of nuclear CRs can still be a gamma-ray 
source with a comparable or even preponderant IC emis- 
sion. Ultimately even SN 1006 might turn out as an ex- 
ample for that case, since its A'h must be quite small . 

4. Conclusions 

From the fact that magnetic field amplification occurs in 
all the young SNRs for which relevant data exist, and given 
the strong theoretical connection between magnetic field 
amplification and efficient acceleration of nuclear CRs, we 
tentatively conclude that the Galactic SNRs are the source 
population of the Galactic CRs. 

Quantitatively, with Cas A, SN 1006, and now Tycho's 
SNR from this paper, there exist three examples where 
nonlinear acceleration theory and X-ray observations of 
filamentary structures give the same morphology and the 
same magnetic field amplification effect: the effective in- 
terior magnetic field energy density is near B'^/^Stt) = 
3.5 X lO^^poVg. The estimates made for the cases of 
Kepler's SNR, RCW 86 and RX J1713.7-3946 confirm 
such a conclusion although less definitely so due to their 
less well-defined parameter values. 

The data suggest that field amplification and thus the 
acceleration of nuclear particles to the highest energies 
that correspond to the knee in the observed CR spectrum 
(Berezhko & Volk 2004b), or even possibly beyond (Bell 
& Lucek 2001; Bell 2004), is a transitory effect in CR 
sources, limited to the early evolutionary phase of SNRs. 
This agrees with theoretical expectations. From acceler- 
ation theory the later phases of SNR evolution, which 
are characterized by decreasing shock speeds, should then 
maintain the lower energy particle population produced 
early on, in a balance between adiabatic expansion and 
continuing acceleration. 
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The highest-energy gamma-rays emitted from the CR 
sources should be observable only during the early phases 
of evolution, roughly at the end of the SNR sweep-up 
phase, before the highest energy particles have already 
left the source. 
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